Proton induced fission of 181 Ta at relativistic energies 
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Total fission cross sections of 181 Ta induced by protons at different relativistic energies have been 
measured at GSI, Darmstadt. The inverse kinematics technique used together with a dedicated set- 
up, made it possible to determine these cross sections with high accuracy. The new data obtained in 
this experiment will contribute to the understanding of the fission process at high excitation energies. 
The results are compared with data from previous experiments and systematics for proton-induced 
fission cross sections. 
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I. INTRODUCTION 



25.85.Ge 



Spallation reactions induced by relativistic protons on 
181 Ta lead to excited target remnants with large fission 
barriers (20-25 MeV) covering a broad range in exci- 
tation energy. The investigation of the fission process 
under these extreme conditions is expected to provide 
relevant information on the dynamics of fission at high 
excitation energies. On the other hand, tantalum and 
tungsten alloys are proposed as optimum materials for 
the construction of spallation neutron sources because 
of their properties under extreme irradiation conditions: 
relatively large neutron production, corrosion resistance, 
and a high melting point. Spallation targets are of inter- 
est in different domains. One of the technologies which 
relies on spallation reactions is that of accelerator-driven 
systems, ADS [l|, which are currently under study as an 
option for nuclear waste incineration. Recently the con- 
struction of the ESS (European Spallation Source) Q 
facility has been approved and different research com- 
munities are awaiting its opening to undertake a wide 
range of experimental programs in material science, bi- 
ology and other scientific disciplines. Tantalum targets 
are also used for the production of exotic nuclei at ISOL- 
type Q facilities and neutrinos [J]. 

Fission may have a significant effect on the perfor- 
mance of a spallation target. Therefore, a good knowl- 
edge of the interaction of protons (commonly used as 
spallation sources drivers) with these materials is manda- 
tory for their characterization. Reactions leading to fis- 
sion are of interest because they contribute to the produc- 
tion of hazardous remnants, in particular gaseous ones, 
such as the isotopes of Kr and Xe. The composition 
of that radioactive inventory, its evolution, the influence 
of these changes to the target performance itself and its 
structural damages can be estimated with state-of-the- 
art models. However, only through an evaluation of nu- 



merical calculations with accurate data it is possible to 
validate these models and improve their reliability for use 
in technical applications. 

Nuclear fission is, indeed, a process which demands a 
complex description of the fissioning nucleus according to 
its excitation energy, angular momentum, fission barrier 
and the time taken to reach the scission point. A first 
approach to the fission description in terms of process 
probability was proposed by Bohr and Wheeler [5[ from 
a purely statistical standpoint. In parallel, Kramers [6| 
introduced a dynamical description of the fission pro- 
cess based on the coupling between internal (excitation 
energy) and collective (deformation) degrees of freedom 
through a dissipation parameter. Based on these ideas, 
Grange et al. [7[ went a step further by including a time- 
dependent solution of the fission width that recently has 
been analytically formulated [8|, Q . Many experimental 
results have already indicated the role of dissipative and 
transient effects in the fission process at high excitation 
energy [TtJ, [H| ■ Manifestations of these effects have also 
been observed in spallation-induced fission reactions 1121 1 , 
in particular with sub-actinide targets 



13] 
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Unfortunately, presently available data related to to- 
tal fission cross sections of 181 Ta above 700 MeV proton- 
beam energy are scarce and show clear discrepancies at 
1000 MeV [14. ll5J. The situation does not improve at 
lower energies, where the available data are more abun- 
dant, but they also present inconsistent results [16h19| . 
in particular between 300 and 500 MeV. Most of these ex- 
periments were performed using passive track detectors 
and only few of them are based on coincident measure- 
ments of both fission fragments [20| . Under such condi- 
tions it seems difficult to unambiguously identify a fission 
channel with a few mb cross section as is expected in this 
case. 

All previous measurements of fission reactions induced 
by protons on 181 Ta were performed using the direct kine- 
matics technique. Therefore, the reaction products have 
a very low kinetic energy, preventing their escape from 
the target. To overcome this difficulty the inverse kine- 
matics technique was utilised in the present work. Sev- 



eral experiments performed at the FRS (FRagment Sep- 
arator) spectrometer using the inverse kinematics tech- 
nique pjl, l2lM25l ] measured the mass and charge of fis- 
sion residues with high precision, in addition to the fis- 
sion cross sections. However, the transmission of the 
fission fragments was limited by the acceptance of the 
spectrometer, and only one fission fragment was mea- 
sured. Therefore, a dedicated experimental setup was 
used in the present work in order to register both fission 
fragments in coincidence with high efficiency and resolu- 
tion m. 

In this work, we show the results of the experiment 
performed at GSI (Helmholtzzentrum fur Schwerionen- 
forschung) aimed at measuring the total fission cross sec- 
tion of 181 Ta induced by protons in the energy range be- 
tween 300 MeV and 1000 MeV taking advantage of the 
inverse kinematics. The dedicated experimental setup, 
enabled the measurement of the cross-sections with high 
precision. Our results will allow for benchmarking the 
state-of-the-art models, and provide data for systematic 
descriptions. 



II. 



EXPERIMENTAL SETUP 



In the present experiment, the 181 Ta nuclei were accel- 
erated using the basic installation of GSI, the UNILAC 
linear accelerator and the SIS-18 synchrotron, up to 300, 
500, 800 and 1000 A MeV with an intensity of the order 
of 10 4 ions/s and a spill duration of 7 s. These beams 
impinged onto a liquid hydrogen target. Due to the kine- 
matics of the reaction, we were able to detect efficiently 
both fission fragments which were emitted in the forward 
direction with large kinetic energies. This capability en- 
abled the use of relatively thick targets, increasing the 
statistics. 



START 



VETO 

I 



DOUBLE SCINT. 



BEAM 



:amB j 



MW 



1 
















1 D 1 




D 




______ — -r- 




1 □ 1 




D 




~~~ — — J 


7T 


SLITS 


MUSIC 1 




MUSIC 2 





FIG. 1: Schematic representation of the experimental setup 
used in the present experiment. 



A sketch of the setup is shown in Fig. [TJ A Multi-wire 
chamber (MW) and thick iron slits were used to collimate 
the 181 Ta beam at the target. A first scintillator detec- 
tor (start) placed upstream of the target, determined the 
beam flux. The target consisted of a liquid hydrogen cell 
(85 mg/cm 2 ) inside a container with 100 /im titanium 
windows and a cryostat for liquefying the hydrogen. The 
target was surrounded by two Multi-sampling Ionization 



Chambers (MUSIC) [27| (200 mm x 80 mm window sur- 
face and 460 mm of active length) which measured the 
energy loss of the tantalum beam particles and that of 
the products of the reaction, respectively. These ioniza- 
tion chambers, having almost 100 % efficiency for the de- 
tection of relativistic heavy nuclei, were used to identify 
reactions of 181 Ta produced in the hydrogen target and 
in any other layer of matter present in the beam line. 
A veto scintillator with a 15 mm diameter hole placed 
just before the target allowed the rejection of beam-halo 
particles and misaligned beam trajectories. 

The two fission fragments were detected independently, 
but in temporal coincidence, by a double paddle scintil- 
lator placed downstream of the target (300 mm x 70 mm 
and 3 mm thickness each paddle). According to the setup 
geometry, two different triggers were used for data acqui- 
sition: The "beam" trigger was provided by the plastic 
scintillator placed upstream of the target in anticoinci- 
dence with the signal of the veto scintillator. 

The "reaction" trigger was produced by the coin- 
cidence between the "beam" trigger and the time- 
coincident signals on both paddles of the double scintil- 
lator placed downstream the target. These two triggers 
provided the measurement of the beam flux together with 
the fission events. The average rates for the "beam" and 
"reaction" triggers were around 10 4 and 700 triggers/s, 
respectively. The "beam" trigger was downscaled to re- 
duce the data acquisition dead time. 



III. 



DATA ANALYSIS 



The identification of the fission events was based on 
the amplitude of the signals recorded by the two MU- 
SIC chambers surrounding the target and the amplitude 
of the signals provided by the two paddles of the dou- 
ble plastic scintillators located downstream of the tar- 
get. With this information, we were able to isolate fis- 
sion events from other reactions channels occurring in the 
hydrogen target. 

In Fig. [2] a scatter plot of the energy losses of ions 
traversing the two MUSICs, before and after the target, 
is depicted using the "reaction" trigger. The events ly- 
ing in the diagonal of this plot correspond to ions which 
kept their atomic number when passing through the tar- 
get. These nuclei, lighter than the primary beam, have 
been produced in nuclear reactions induced by 181 Ta 
projectiles in any layer of matter situated upstream of 
the hydrogen target. The dominant 181 Ta spot of non- 
interacting beam particles is clearly visible at the top, 
near channel 3800 on the vertical axis. The vertical group 
below the beam spot corresponds to residual fragments 
produced in the interaction of 181 Ta with hydrogen. In 
this group, events inducing high and small energy loss 
signals, correspond to residual heavy nuclei and emitted 
light nuclei from evaporation processes. 

Since the energy loss of nuclei is proportional to their 
atomic number squared (Z 2 ), fission fragments are ex- 
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FIG. 2: Scatter plot of the amplitudes of the signals registered 
with the two MUSIC detectors. Nuclei lighter than 181 Ta pro- ^ 

duced in reactions before MUSIC 1 appear in the diagonal 
region. In the vertical line the 181 Ta spot corresponding to 
non-interacting beam particles is represented, and below the 
events corresponding to reactions in the target. The box en- 
closes the fission region and the color code represents counts 
on the logarithmic scale. 
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pected to produce energy loss signals corresponding to 
about half of the value obtained for the primary beam 
(AE f j. ex Z\ +Z$ = Zl eam /2). Therefore, fission prod- 
ucts should be located around channel 1800 on the MU- 
SIC 2 energy loss axis. To count for the fission events, 
nfi ss , a condition in the scatter plot shown in Fig. [2] was 
applied selecting the region where the fission products 
are expected. Focusing on the selected region, indicated 
by the rectangular area in Fig. [3J the fission events were 
identified combining the amplitude (energy loss) of the 
signals recorded by the two paddles of the double plastic 
scintillator. 

In Fig. [3] the amplitudes of the signals registered by 
both plastic scintillators in temporal coincidence (at 1 A 
GcV in upper panel and 300 A MeV lower panel), us- 
ing the "reaction" trigger, are represented in a scatter 
plot. Due to the charge splitting of the fission process, 
fission events are expected to populate the diagonal band 
in this figure, and are separated from other much more 
abundant reaction channels. This fission region only rep- 
resents a small part of the plot statistics since the fission 
probability is rather small. For this reason, fission events 
could only be properly identified by a detection setup 
enabling the identification of different reaction channels. 

In order to provide an accurate measurement of the fis- 
sion cross section, we evaluated the background which re- 
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FIG. 3: Scatter plot of the energy loss signals provided by 
the two paddles of the double plastic scintillator placed down- 
stream the target with the reaction trigger (Upper panel 1 A 
GeV Full target - Lower panel 300 A MeV Full target). Both 
plots are normalized to the same number of counts to high- 
light the influence of different reaction channels. The color 
code represents counts on the logarithmic scale. 



mains in the fission region due to simultaneous break-up 
and evaporation processes. To evaluate this background 
we used Fig. 01 where the energy loss provided by the two 
paddles for the double plastic scintillator at 1000 A MeV 



is represented by selecting only events compatible with a 
fission signal in the MUSIC detectors (rectangular area 
in Fig.|2). 

In this figure, intermediate mass fragments (IMFs) 
produced in simultaneous break-up reactions may pop- 
ulate the fission region (dotted contour in Fig. HJ) . The 
evaluation of this break-up background was performed 
via dividing the fission region into slices as shown by the 
thin rectangles in Fig. [4] Each slice was then projected 
along its longitudinal dimension (insets in Fig. |4|) , which 
clearly enhanced the profile of the contributions com- 
ing from background (left peak) and fission (right peak) . 
Gaussian fits to each of the two contributions defined the 
correction for the break-up background suppression. On 
the other hand, evaporation residues could also populate 
the edges of the fission region along an axis AE\ +AE2 
(dashed line in Fig. 2]) defined by the sum of the signals 
of the two scintillators. To overcome this problem, the 
region profile (dotted contour in Fig. \Q was projected 
onto this AEi +AE2 axis to evaluate this contribution 
by means of gaussian fits as shown in Fig. [5] 

The number of measured fission events nfi SS corre- 
sponds then to the number of events in the fission region 
corrected by the overlapped simultaneous break-up and 
evaporation background contributions. 
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FIG. 5: Projection of the fission region defined in Fig. [4] on 
the AEi +AE2 axis represented in the same figure. The con- 
tribution of evaporation residues was evaluated by means of a 
gaussian fit (dotted line) and subtracted from the total con- 
tribution (dashed line). 
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FIG. 4: Same as the upper panel in Fig.[3]but conditioned by 
the fission selection from Fig [5] The different contours and 
insets illustrate the background suppression method used to 
identify fission events as explained in the text. 

Fission yields (Yfi SS ) were obtained from fission event 
measurements corrected by the background (n/, ss ) and 
additional effects such as the secondary reactions of the 
fragments in the target (fd) and the geometrical accep- 



tance of the experimental setup {J ge o)i according to the 
following equation: 



Y, 



fiss 



Tlfiss Jd ' Jg 



(i) 



Secondary reactions of the fission fragments in the 
target were evaluated using the Karol's microscopic 
model [28| and amounted to less than 2.5 % for full tar- 
get and less than 0.5 % for empty target measurements 
at 1000 A MeV. Geometrical constraints were also con- 
sidered to evaluate the efficiency of the detection setup. 
Fission products emitted close to the double paddle scin- 
tillator gap had a probability of passing through it or 
through the same paddle. A Monte Carlo simulation 
based on the liquid-drop-model calculation of the post- 
scission kinetic energy of the fission fragments [29[ was 
performed to evaluate the ratio of fission product losses 
due to the geometry of the setup. Taking into account the 
dispersion of the beam as measured with the Multi-wire 
chamber detector and the alignment, the coulomb force 
between both fragments and the distance from the centre 
of the hydrogen target to the double plastic scintillator, 
we calculated the perpendicular dimensions of the fission 
fragments distribution in the double plastic scintillator 
detection plane. The acceptance of these scintillators 
(300 mm x 140 mm and 1 mm gap) allowed us to evaluate 
the losses. Thus, the resulting yield was corrected by a 
geometrical factor (f geo ) having a value smaller than 10 % 
at 1000 A MeV and decreasing with the beam energy. 



To determine the number of projectiles (rib), we used 
the first MUSIC to identify tantalum among other nuclei 
that have been created in other layers of matter placed 
in the beam line before the target, as shown in Fig. [5] 
The sum of the Z=73 ions identified according to this 
procedure using the "beam" and the "reaction" trigger 
corrected by the downscaling factor provided the total 
number of projectiles. 
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sources of systematic uncertainty for the different cor- 
rection factors were also evaluated. The systematic un- 
certainty of the geometrical correction factor (f ge o) was 
estimated to be smaller than 5 %. The evaluation was 
done by changing the size of the double plastic scintil- 
lator gap and the beam profile in our simulation. The 
value of the systematic uncertainty of the correction fac- 
tors due to the beam attenuation (J a ) and the secondary 
reactions of the fission fragments (fd) were smaller than 
1 % and almost the same for all energies. Due to the 
relatively large number of recorded fission events the sta- 
tistical uncertainties were below 1.5 %. Statistical and 
systematic uncertainties, others than the ones associated 
to the beam intensity and target thickness, for the mea- 
surements with the full and empty target are presented 
in Table U and Table HU 



TABLE I: Statistical (e 3ta t.) and systematic uncertainties due 
to the identification of fission fragments (e(n/i ss )), geomet- 
rical acceptance (e(f geo )) and the attenuation (e(f a )) of the 
beam affecting our measurements. 



Energy 



^stat. 



f(n /lss ) e(fge ) e(f a ) e(£j) 



1000 A MeV 


0.38% 


6.01% 


4.32% 


0.43% 


0.22% 


800 A MeV 


0.39% 


7.16% 


3.72% 


0.41% 


0.22% 


500 A MeV 


0.45% 


7.84% 


2.80% 


0.39% 


0.22% 


300 A MeV 


0.33% 


9.51% 


2.05% 


0.38% 


0.22% 



FIG. 6: Energy Loss signals provided by the MUSIC 1. The 
region between the dotted red lines corresponds to 181 Ta. 

Due to the relatively large thickness of the target, the 
number of projectiles was corrected by a factor f a (< 5 % 
for full target and < 1 % for empty target at 1000 A 
MeV) taking into account the attenuation of the beam 
intensity along the target. The value of this factor was 
also evaluated using Karol's model. 

To correct for reactions taking place in the target win- 
dows (namely Ti) which surrounded the liquid hydrogen, 
fission yields determined with the empty target following 
the same analysis procedure were subtracted from the 
fission yield obtained with the full target. Finally, the 
respective fission yields were normalized to the number 
of projectiles and the number of nuclei in the target per 
unit area (N t ) to determine the total fission cross section 
according to the following expression (with JVj = nyfa): 



Y. 



full 



fiss 



N, 



full 



yempty ' 

fiss 
jy-empty 



1 

N't' 



(2) 



Particular attention was paid to the evaluation of the 
corresponding uncertainties. The main sources of sys- 
tematic uncertainties were the identification of fission 
events ranging from 2 to 10% (e(n/j SS )), the beam in- 
tensity (m 5%) and the target thickness (« 4%). The 



TABLE II: Same as Table U but for empty target measure- 
ments. (*At 500 A MeV the independent number of fission 
could not be determined but the yield normalized to the num- 
ber of projectiles.) 



Energy 



£stat. 



sjnfias) e(f geo ) e(fq) e(f d ) 



1000 A MeV 


1.45% 


1.95% 


4.32% 


0.05% 


0.04% 


800 A MeV 


1.05% 


4.56% 


3.72% 


0.05% 


0.04% 


500 A MeV* 


- 


6.71% 


2.80% 


0.05% 


0.04% 



300 A MeV 0.96% 7.08% 2.05% 0.05% 0.04% 



IV. RESULTS 

Using the method described in the previous sections, 
we have measured with high precision the total fission 
cross section of 181 Ta induced by protons at 300, 500, 800 
and 1000 A MeV. The results obtained for each energy 
are presented in Table Mil The magnitude of the mea- 
sured cross sections is rather small and strongly decreases 
for the lower beam energies. The associated uncertain- 
ties are also rather small (« 10 %) but increase for the 
lowest energies (« 18 %) since the smaller fission cross 
sections complicates the identification of fission events. 



TABLE III: Total fission cross sections determined in this 
work 



Energy 


Fiss. cross section 


Stat, uncert. 


Syst. uncert. 


E 


(A MeV) 


(mb) 


(%) 


(%) 


c 


1000 


20.17 ± 2.19 


0.46 


10.85 


*■» 


800 


13.09 ± 1.62 


0.32 


12.34 


o 

0) 


500 


7.53 ± 1.40 


0.51 


18.54 


to 


300 


6.55 ± 1.00 


0.48 


15.21 


in 
o 



In Fig. we present the cross sections obtained in 
this work as solid points compared to previous measure- 
ments by different authors. In this figure, we also present 
predictions obtained from the systematics established by 
Prokofiev some years ago [30| (dashed line). 

From the analysis of the previously measured cross 
sections, one can identify some clear discrepancies. 
At the highest energies, one can find two rather dis- 
crepant measurements around 670 MeV by Konshin et 
al. [H (14.0 ± 1.9 mb) and by Baranovskiy et al. [HI 
(8.0 ± 2.5 mb). At 800 MeV it exists a single mea- 
surement by Yurevich et al. 14 1 and at 1000 MeV one 



finds again two discrepant values obtained by Yure- 
vich et al. (15.65 ± 5.4 mb) and Bochagov et al. [l5[ 
(27.0 ± 1.5 mb). Our results are in very good agreement 
with the measurement of Yurevich et al. at 800 MeV 
and within the error bars at 1000 MeV, solving the exist- 
ing discrepancy in this energy range. Moreover, we also 
confirm the predictions estimated by the systematics of 
Prokofiev. 

In the energy range between 300 and 600 MeV, we 
can also observe important discrepancies between differ- 
ent measurements. Around 300 MeV the data obtained 
by Yurevich et al. (5.2 ± 1.6 mb) and Konshin et al. 
(2.6 ± 0.4 mb) differ by a factor two. Around 400 MeV, 
the measurements by Yurevich et al. (5.79 ± 1.78 mb) 
and Konshin et al. (4.7 ± 0.7 mb) are in rather good 
agreement. However, the measurement by Konshin et al. 
is significantly smaller than the one obtained from the 
systematics of Prokofiev (7.60 mb). The measurement by 
Yurevich et al. could be compatible with the systematics 
owing to its large uncertainty. Finally, around 500 MeV 
the measurements by Konshin et al. (8.3 ± 1.1 mb) and 
Yurevich et al. (5.59 ± 1.72 mb) also differ by a large 
factor. Our measurements at 300 and 500 MeV are con- 
sistent with the estimated values from the Prokofiev for- 
mula, and, confirm the largest values of the cross sections 
measured in this region. 

From this analysis we can conclude that our data con- 
firm the measurements by Yurevich et al. and the cross 
sections estimated by the systematics above 700 MeV. At 
lower energies our measurements clarify the discrepancies 
existing until now. In the energy range between 300 and 
600 MeV, our data favor those measurements presenting 
the highest cross sections. Moreover our data confirm 
the predictions obtained by the systematics of Prokofiev 
over the entire energy range covered by this work. 
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FIG. 7: Fission cross sections measured in this work (solid 
circles) in comparison to previously measured data and esti- 
mates obtained from the systematics established by Prokofiev 
(dashed line). 



SUMMARY AND CONCLUSIONS 



In this work we have investigated the proton induced 
fission of 181 Ta in inverse kinematics at 300, 500, 800 
and 1000 A MeV. The combination of the inverse kine- 
matics technique with a highly efficient detection setup 
made it possible to determine the total fission cross sec- 
tions with high accuracy. The coincident measurement 
between both fission fragments and their identification 
from the rough determination of their atomic number al- 
lowed to clearly identify and separate the fission events 
from other reaction channels. This selection is shown to 
be extremely useful at lower energy when the fission cross 
section is small and complicates the identification of this 
reaction channel. 

The new data have completed the scarce number of 
measurements existing above 700 MeV. At intermediate 
energies, the quality of the new data enabled to clarify 
previous results. Moreover, these new data are overall 
in good agreement with the systematics established by 
Prokofiev over the entire energy range. The results will 
also allow to benchmark different state-of-the-art models 
describing fission and transport simulation codes used in 
the design of spallation targets. 
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